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Abstract 

The recently published cosmological bound on the absolute neutrino masses obtained from the 
Wilkinson Microwave Anisotropy Probe (WMAP) data has important consequences for neutrino 
experiments and models. Taken at face value, the new bound excludes the determination of the 
absolute neutrino mass in the KATRIN experiment and disfavors a neutrino oscillation interpre- 
tation of the LSND experiment. Combined with the KamLAND and Super-K data, the WMAP 
bound defines an accessible range for the neutrinoless double beta decay amplitude. The bound 
also impacts the Z-burst annihilation mechanism for resonant generation of extreme-energy cosmic 
rays on the cosmic neutrino background in two ways: it constrains the local over-density of neutrino 
dark matter which is not helpful, but it also limits the resonant energy to a favorable range. In 
R-parity violating SUSY models, neutrino masses are generated by trilinear and bilinear lepton 
number violating couplings. The WMAP result improves the constraints on these couplings by an 
order of magnitude. 
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I Introduction 



With the recently pubhshed first data of the Wilkinson Microwave Anisotropy Probe (WMAP) [Ij on 
the cosmic microwave background (CMB) anisotropics the age of precision cosmology has arrived. A 
flat, vacuum-energy dominated cold dark matter (ACDM) universe seeded by nearly scale-invariant 
Gaussian primordial fluctuations appears to be firmly established as the standard cosmology. More- 
over, when combined with additional CMB data-sets (CBI, ACBAR) j2j and observations of large 
scale structure from the 2dF Galaxy Redshift Survey (2dFGRS) jS] to lift degeneracies, the WMAP 
data offers the potential of testing various extensions and sub-dominant components in the ACDM 
model, such as small non-flatness, quintessence, possible tensor-gravitational wave modes, and a mas- 
sive cosmic neutrino background (Ci^B). Investigation of the latter has most important consequences 
for terrestrial physics experiments exploring the neutrino sector. 

The power spectrum of early-Universe density perturbations is processed by gravitational instabilities. 
However, the free-streaming relativistic neutrinos suppress the growth of fluctuations on scales below 
the horizon (approximately the Hubble size c/H{z)) until they become nonrelativistic at 2; mj/3TQ ~ 
1000 {nij /eV) . When the amplitude of fluctuations is normalized to the WMAP data, the amplitude 
of fluctuations in the 2dFGRS places significant limits on the contribution of neutrinos to the energy 
density of the universe, 

^-h^ = I^fl < 0-0076 (95%C.L.), (1) 
9o.5 eV 

which translates into 

^mi< 0.71 eV (95%C.L.). (2) 

i 

The new mass bound ((SJ impacts most directly four-neutrino mass models constructed to accommodate 
the LSND evidence for oscillation. Such models require the heaviest neutrino mass to be ~ 1 eV, and 
so at face value are disfavored by the new result [U IHj. However, there are several loopholes in 
the argument against an ~ 1 eV sterile neutrino. If there is only one isolated "heavy" sterile as in 
the 3+1 model, then the WMAP/2dF data at face value allow the Aml^^Y) region up to 0.5 eV^, 
whereas relaxing the WMAP/2dF bound from 0.71 eV to 1 eV allows virtually the entire LSND 
region to co-exist. In a 2+2 model, there are two heavy mass eigenstates, and the WMAP/2dF data 
at face value limit AmLSND *° 0.1 eV^. Still another possibility, not yet explored to the best of our 
knowledge, might be to model the heavier neutrinos as decaying to light flavors plus a light boson, 
with a lifetime much less than the age of the Universe at structure formation. In such a model, the 
decay products would be free-streaming particles with masses well below the WMAP bound. Relevant 
to this discussion is the limit from Big Bang Nucleosynthesis (BBN) that neutrinos beyond the 

three active could not have been in thermal equilibrium already at the BBN temperature ~MeV, long 
before the epoch of structure formation. So the more serious constraint for the sterile neutrino is the 
BBN limit. Overcoming this BBN limit automatically immunizes the sterile against the WMAP/2dF 
bound IS], since the depopulated states at BBN are not populated at a later time. One way to evade 
thermalization at the BBN epoch is via a tiny lepton asymmetry [H]. There are several other ways, 
conveniently summarized in [S]. In conclusion, MiniBooNE is still required to settle the fate of the 
sterile neutrino jlOj . 

From here on we focus on the consequences of the new WMAP bound for three-neutrino models. 
It was previously noted that there are potentially four independent approaches for measuring 
the absolute neutrino mass. These are large-scale structure studies measuring the total mass in the 
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Figure 1: Implications of the WMAP neutrino mass bound for the mass of the heaviest neutrino 7713. Here we 
take the best-fit value for Arnl^^ = 3 x 10^^ eV^; Airi^^^ is too small to be relevant. 



Cz^B (as reported by WMAP), the Z-burst method measuring individual masses in the Ci^B, and the 
terrestrial measurements of the tritium end point spectrum and neutrinoless double beta decay rate. 
Of course, the results of these approaches are correlated in the sense that the experiments all attempt 
to determine the same neutrino masses. We will examine the impact of the new WMAP bound on 
the future of the other three approaches. 

Neutrino oscillation studies have established three important facts of relevance here. The first is that 
the two mass-squared differences are small compared to the WMAP limit. Thus, when the WMAP 
limit is saturated, the three neutrinos are nearly degenerate in mass, and we have 

mi<0.24eV (95%C.L.) (3) 

for each neutrino mass. The second important fact is that oscillation studies provide a lower bound 
on the heaviest neutrino mass, given by the minimum y'Am^^ ~ 0.03 eV. Thus, we may write 

0.03 eV < ma < 0.24 eV (95% C.L.) , (4) 

which shows the remarkable fact that knowledge of the heaviest neutrino mass (which we shall always 
denote by m^) is now known to an order of magnitude! A plot of the total neutrino mass versus 1713 
is shown in Fig. ^ The relation is linear, "^^rrii = 3 771.3, except very near the smallest allowed ms, 
~ \/ Am^^-jjj. The third important fact is that the three angles parameterizing the unitary flavor-mass 
mixing- matrix, Uai, are well known. The one CP-violating Dirac phase and two CP- violating Majorana 
phases are not known. The angles and phases will be important when we look at neutrinoless double 
beta decay. 

Absolute neutrino mass bounds also constrain all entries in the neutrino mass matrix in flavor space 
due to unitarity. This results in bounds on couplings in theories with lepton number violation As 
an example, we derive bounds on parameters of the R-parity violating {B/p ) SUSY model, improving 
them by one order of magnitude over the existing values. 
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a) 



b) 



c) 



Figure 2: Neutrino mass spectra for the three neutrino case: a) degenerate, b) hierarchical and c) inverse 
hierarchical spectrum. 



II Tritium beta decay 

The mass to be inferred from /3-decay is = Y,j Pej?' m'^j- The KATRIN project JJj plans to start 
taking data in 2007. The sensitivity aim after three years of measurement is 0.08 eV^ at 1 o" accuracy. 
This may be improved to 0.05-0.06 eV^, when optimizing the data point distribution and resolution, 
which implies a final sensitivity of m,^^ to be 0.4 eV at 3 a. Thus, the reach of this experiment includes 
only the nearly mass-degenerate neutrino case, for which unitarity allows one to write m^^ = m^. 

Comparing the KATRIN reach to the WMAP limit in Eq. © , one comes to the unfortunate conclusion 
that a positive signal is unlikely. 



Ill Neutrinoless double beta decay 

The mass inferred in neutrinoless double-/? decay is 

niee = \^Ulimi\ . (5) 

i 

Here one needs the neutrino mixing parameters explicitly. The most recent analysis of atmospheric 
neutrino data yields 

10^^ eV^ < Am^^^ < 5 • 10~^ eV^ (6) 

and 

sin^ 20atm > 0.8. (7) 

On the other hand, a recent evaluation of solar neutrino data including the KamLAND reactor exper- 
iment TS' inferred 

5 • 10"^ eV^ < Am^^^ < 1.1 ■ 10"^ eV^ (8) 
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and 

0.3 < tan^ 6'sun < 0.8 . (9) 

Thus, the LMA solar solution is confirmed. The neutrino mixing matrix is seen to be "bi-large." It is 
also known that |C/e3p ~ 0, which means that the third mixing angle is negligibly small jl6j . 

The cases of degenerate, hierarchical, and inverse hierarchical neutrinos (see Fig. ^ must be considered 
separately (for a detailed discussion, see e.g. ^ZJ). The WMAP limit is sufficiently large that it impacts 
only the case of degenerate neutrinos. 



• Degenerate Neutrinos: mi ~ m2 — rris. With |C/e3p ~ 0, one has a mass proportional to 
\Uei + U^2\ = I cos2 9 + e^*^ sin^ 9\, which, upon extremizing the unknown phase, leads to 

cos 26lsun"T'3 < "Ige < "13 . (10) 

Inputting the new WMAP bound, and the solar angle, one gets 

0.1 ms < mee < 0.24 eV. (11) 

• Hierarchical neutrinos: mi <^ m2 <C m^ and Am"^^^ = Amfg • Here a lower limit is obtained by 
taking |C/e3p = and mi = 0. The result is 

mee > VAm2^„sin2 e^nn = 2 • 10-^ eV, (12) 

and mee ^ "^3 ~ \/^"Sm 

^ 0.07 eV. 

• Inverse hierarchical neutrinos: mi <C m2 — m^ and Am"^^^ = Am2-^. The situation is analogous 
to the degenerate case, but with the scale of 7713 fixed by the atmospheric neutrino evidence, 
rather than the WMAP result. One gets 

cos 2e,^r,\J ^rnl^^ < mee < \J Am2^^ , (13) 

i.e. 

3 • 10^3 eV < mee < 0.07 eV . (14) 



In summary, neglecting unnatural cancellations due to a conspiracy of 5, mi and mixing angles, the 
predicted range of mee is given by 

2 • 10"^ eV < mee < 0.24 eV. (15) 

Fortunately, the whole region can be covered by the most ambitious double beta decay proposals ^Hl 
(for an overview of the experimental status see [TO])- The lower limit is not impacted by the WMAP 
result, whereas the upper limit comes directly from the WMAP data. The central value of the recent 
discovery claim of the Heidelberg-Moscow experiment [201) "^ee = 0.391q'3| eV, exceeds the WMAP 
bound, but the reported lower range does not (this fact has been pointed out already in ref. _H]). 
We point out, though, that double beta decay mechanisms other than the standard neutrino mass 
mechanism are not affected by this bound. A particular interesting possibility to accommodate the 
Heidelberg-Moscow result involves singlet neutrinos propagating in large extra dimensions in which 
case a mechanism decorrelating the neutrino mass eigenstates from the double beta decay amplitude 
is operative |2j|. Exchange of superpartners in R-parity violating SUSY, leptoquarks, or right-handed 
W bosons constitute other possibilities to account for a sizable neutrinoless double beta decay signal 
(for a review see |22j). 
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Figure 3: Resonant energies for different neutrino mass eigenstates in the Z-burst model as a function of the 
largest neutrino eigenmass m^. The Arris's used here are the same as in Fig. 1. 



IV The Z-burst model for EECR's 



The Z-burst mechanism j^Hj generates extreme-energy cosmic rays (EECR's) by resonant annihilation 
of a EECR neutrino on the Cz^B neutrinos. The resonant energy is 

The decay products of Z-bursts include on average two nucleons and, from ten neutral pions, twenty 
photons. The decay multiplicity is ~ 30. The nucleons lose / ~ 20% of their energy for each 
A ~ 6 Mpc traveled in the CMB, so the average energy of a secondary nucleon arriving at Earth from 
distance D is 

lO^ieVx (0.8)^/6 ^^P^ 
~ (m./0.1eV) • 

The photons have shorter absorption lengths, except above 10^^ eV, and so are not expected to con- 
tribute much. For a neutrino mass in the range of Eq. @, the mechanism is optimized: a larger mass 
would move Z-burst secondaries down below the GZK energy ~ few x 10^^ eV where the "background" 
of normal EECR events is huge, whereas a smaller mass would move the resonant energy beyond the 
reach of any realistic neutrino flux. The Z-burst resonant energies as a function of the heaviest neu- 
trino mass 7713 are shown in Fig. |31 Note that over most of the allowed range, all three neutrinos 
contribute to annihilation with a resonant energy within a factor of two of each other. 

In the simplest approximation, the spectrum of arriving nucleons is 

dN I dN dD ^ -I 

r -X X oc (18) 

dE L>2 dD dE ^ ' 



from sources uniformly distributed out to 



^"^■^-^ Ml-/) ■ (1^' 
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with a pileup at Eqzk resulting from all primaries originating beyond this distance. The 1/E spectrum 
extends from Eqzk out to the maximum nucleon energy ~ E^/m ~ 1021 (0^) eV. More realistic 
simulations including energy-loss processes, cosmic expansion, and boosted Z-boson fragmentation 
functions give a softer spectrum, but a characteristic feature of the Z-burst mechanism remains that 
the super-GZK spectrum is considerably harder than the sub-GZK spectrum having power-law index 
-2.7. 

What is not known is whether Nature has provided the large neutrino flux at E^ to allow an appreciable 
event rate in future EECR detectors. It is conceivable, although unlikely, that the flux is so large that 
present EECR events are initiated by Z-bursts. A recent analysis 24 of this possibility gave a best 
fit with nil, = 0.261q'2^ eV, nicely consistent with the WMAP bound. Another analysis |25j fits the 
EECR spectrum down to the ankle with Z-burst generated events and a neutrino mass of 0.07 eV, 
again in accord with the WMAP bound. The flux requirements for the Z-burst mechanism can be 
ameliorated if there is an over-density of relic neutrinos, as would happen if (i) there was a significant 
chemical potential, or (ii) neutrinos were massive enough to cluster in "local" structures such as the 
Galactic SuperCluster. Large chemical potentials have been ruled out recently |2H], and this exclusion 
is confirmed by the WMAP data. Local clustering has been studied [27j, with the conclusion being 
that a significant over density on the SuperCluster scale requires a neutrino mass in excess of 0.3 eV. 
Such a mass is marginally allowed by the new WMAP/2dF limit. 



V WMAP neutrino mass bound on Ijtp SUSY 



Supersymmetry without i?-parity [2Sj provides an elegant mechanism for generating neutrino (Majo- 
rana) masses and mixings. In these models, there are mainly two sources of neutrino mass generation. 
In one scenario, products of trilinear A and/or A' couplings generate a complete neutrino mass matrix 
through one-loop self-energy graphs |29| In the other scenario, the bilinear i?-parity- violating 

terms induce sneutrino vacuum expectation values (VEV's) allowing neutrinos to mix with the neu- 
tralinos. 

The L- violating part of the superpotential can be written as 

y^RPV = ^KjkLiLjE^ + X[ji^LiQjDl + fiiLiHu, (20) 

where i, j and k are quark and lepton generation indices. In Eq. (PU)) . Li and Qi denote SU(2)-doublet 
lepton and quark superfields, E^ and are SU(2)-singlet charged lepton and down-quark superfields, 
and Hu is the Higgs superfield responsible for the mass generation of the up-type quarks, respectively. 
There are 9 A-type (due to an antisymmetry in the first two generation indices), 27 A'-type and 3 Hi 
couplings. Stringent upper limits exist on all these couplings from different experiments [311 132j . 

We first consider the effects of the A' interactions. The relevant part of the Lagrangian can be written 
as 



^X' — Kjk 



dkPhi'idjL + lyfPLdjdlji 



+ h.c. (21) 



Pl is the left-helicity projector. Majorana mass terms for the left-handed neutrinos, given by Cm = 
— ^fnu^^ii^Lii^jii' + h.c, are generated at one loop. Figs. |1] show the corresponding diagrams. The 
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masses so induced are given by 



where fix) = {xlnx — x + l)/{x — 1)^. Here, m^,; is the down quark mass of the i*^ generation 
inside the loop, is some kind of an average of dLi and dju squark masses, and A^c = 3 is the 
color factor. In deriving Eq. (|22() . we assumed that the left-right squark mixing terms in the soft 
part of the Lagrangian are diagonal in their physical basis and are proportional to the corresponding 
quark masses, i.e. Amf^^{i) = rrid-m^ . The small effect of quark mixing is neglected in order not to 
complicate the discussion unnecessarily. 

With A-type interactions, one obtains exactly similar results as in Eqs. (|2H) and (|22|) . The quarks 
and squarks in these equations will be replaced by the leptons and sleptons of the corresponding 
generations. The color factor Nc = 3 would be replaced by 1. We do not explicitly write them down. 

For numerical purpose, we have assumed the mass of whatever scalar is relevant to be 100 GeV 
throughout, to be consistent with common practice and, in particular, to compare with the old bounds. 
While for sleptons this sounds a reasonable approximation, for squarks the present lower limit, even 
in I^p scenarios, is around 250 GeV [HSI- In any case, for different squark masses one can easily derive 
the appropriate bounds by straightforward scaling. It should be noted that the product couplings 
under consideration contribute to charged lepton masses as well, but with the present limits those 
contributions are too small to be of any relevance. The resulting bounds are 

A^33A^,33 < 3.6 • 10-^ A^32A^,23 < 8-9 • 10-^ A^22A^'22 < 2.2 • 10-^ 
Ai33Ai33 < 6.3 • 10-^ A,32Ai23 < 1.1 • 10-^ A,22Ai22 < 1.7 • 10-^ (23) 

There is one combination which receives a more stringent limit from conversion in nuclei j34j . 
namely A']^22A222 < ^-"^ ' lO"''- The chirality flips in Figs. explain why with heavier fermions inside 
the loop the bounds are tighter. For this reason, we have presented the bounds only for j, k = 2,3. 

Next we consider the bilinear terms. Such terms lead only to one massive eigenstate as a result 
of tree level mixing between neutrinos and neutralinos. The induced neutrino mass |35j is given by 
m ~ Assuming the Higgsino mixing parameter /i = 100 GeV, one obtains 

m/fi < 1.5 • 10"^ (24) 

The bounds in Eqs. ()23|) and (|24() obtained using the recent WMAP bound are more stringent than the 
existing ones by one order of magnitude, precisely to the extent that the WMAP data have improved 
the absolute neutrino mass bound. 

We make a note in passing that even our improved bounds on trilinear couplings do not invalidate 
the F^p SUSY search strategies proposed by the authors of ,36;,. Their suggestion is that at the 
Tevatron collider the production and decay of sparticles would occur in R-parity conserving modes 
except that the neutralino LSP would decay via I^p channel into multi-5 and missing energy final 
states constituting the signal. 
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(22) 
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VI Conclusions 



We have discussed implications of the WMAP neutrino bound on future neutrino mass studies, in- 
cluding Tritium beta decay, neutrinoless double beta-decay, and the Z-burst mechanism for EECR's. 
We have shown that the Tritium beta decay project KATRIN is unlikely to measure an absolute 
neutrino mass, and that the WMAP bound in combination with the neutrino oscillation data defines 
a predicted range for the double beta-decay observable rriee, which is accessible in the most ambitious 
proposed experiments. The WMAP bound also impacts the Z-burst mechanism for cosmic rays above 
the GZK cutoff. It constrains local over-densities, but it also limits the resonant energy to a favorable 
range. 

Turning to model building, WMAP constrains theories with AL = 1 lepton number violation, since 
in these theories Majorana neutrino masses are generated radiatively. Taking B/p SUSY as our 
example, we have derived the upper limits on many individual and product couplings of the A- and 
A'-types, and also the bilinear /ij terms, from their contribution to neutrino masses. Using the recent 
WMAP bound the limits have been improved by an order of magnitude. Finally, we remark that the 
new WMAP bound on neutrino mass coincides nicely with the one arising from the requirement of 
successful baryogenesis in the context of the neutrino see-saw model |37) . 



Note Added 



The new mass upper-limit expressed in Eqs. @, ©, and (0J) depends on priors. In particular, the 
bound depends on the inclusion of Lya data to estimate the power spectrum of intervening hydrogen 
clouds. Given the complexity of the Lya analysis, some have questioned the reliability of this data 
set. Without the Lya prior, the WMAP mass limit is relaxed to 0(1) eV [^H]. Accordingly, the upper 
dot-dash lines in Figs. Hand 01 and the upper bounds in Eqs. Ijllj) . (|^ . and would also relax 
by ~ 40%. 
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Figure 4: The A'-induced one loop diagrams contributing to Majorana masses for the neutrinos. The A-induced 
diagrams are analogue with sleptons propagating in the loop. 
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